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We calculate effects of an applied helically symmetric potential on the low energy electronic 
spectrum of a carbon nanotube in the continuum approximation. The spectrum depends on the 
strength of this potential and on a dimensionless geometrical parameter, P, which is the ratio of 
the circumference of the nanotube to the pitch of the helix. We find that the minimum band gap of 
a semiconducting nanotube is reduced by an arbitrarily weak helical potential, and for a given field 
strength there is an optimal P which produces the biggest change in the band gap. For metallic 
nanotubes the Fermi velocity is reduced by this potential and for strong fields two small gaps appear 
at the Fermi surface in addition to the gapless Dirac point. A simple model is developed to estimate 
the magnitude of the field strength and its effect on DNA-CNT complexes in an aqueous solution. 
We find that under typical experimental conditions the predicted effects of a helical potential are 
likely to be small and we discuss several methods for increasing the size of these effects. 



PACS numbers: 73.22.Dj, 73.43. Cd, 73.63.Fg 
I. INTRODUCTION 

In recent years it has become a common practice to 
functionalize CNTs with water soluble, high molecular 
weight, quasi-linear molecules such as synthetic polymers 
and DNA.^ ^ These molecules bind to individual CNTs 
via van der Waals forces and form robust and stable 
complexes with the nanotube.'^ The resulting complex 
is easily dispersed in an aqueous solution because the hy- 
drophobic nanotube is screened from the water by the 
wrapping molecule, while the hydrophilic regions of the 
DNA or polymer are free to interact with the solvent. 
Forming such complexes between CNTs and quasi-linear 
molecules is the only known way to disperse nanotubes 
in aqueous media without using surfactants or chemi- 
cally modifying the CNT. Several practical applications 
have been realized as a result of this technique. DNA 
wrapped CNTs can be sorted by diameter using density 
gradient ultra-centrifugation, and the resulting mixture 
filtered to obtain a solution comprised almost entirely 
of one nanotube species.i^ DNA wrapped CNTs are es- 
pecially well suited for biological applications, such as 
cellular markers,^ which are not possible with surfactant 
dispersed CNTs because surfactant molecules generally 
destroy biological systems. 

Simulations and experimental observations indicate 
that polymers and DNA can wrap a round th e exterior of 
a CNT in an ordered, helical fashion t^jg paper 
we will study the effect of a helical potential on the sin- 
gle particle energy spectrum of a CNT. From a practical 
perspective, such a calculation is needed because of the 
large number of experiments performed on DNA-CNT 
complexes and changes in the band structure induced by 
the helical potential must be taken into account when 
interpreting experimental data. Also, while many poly- 
mers easily dissociate from the nanotube upon a change 
in solvent,^ DNA is much more difficult to remove.'^ For 
device applications, it will be useful to know if a DNA- 
CNT complex can be substituted for a pristine CNT 



without a change in response or loss of function. From a 
more fundamental perspective, our work provides scaling 
relationships relating changes in the band structure to 
structural parameters such as the nanotube radius and 
the pitch of the helix. 

Two related studies have recently been pubhshed. In 
Ref. (fTTj) . the authors modeled the nanotube as a free 
electron gas confined to the surface of a cylinder and 
governed by the non-relativistic Schrodinger equation. A 
helical potential was introduced as a series of delta func- 
tions, similar to the Kronig-Penney model. This study 
concluded that the total electronic energy is an oscil- 
latory function of the pitch, with several local minima 
indicating preferred wrapping angles. However, it has 
been shown that the low energy electronic structure of 
a CNT is better described by the D irac Hamiltonian for 
a massless relativistic particle .^i^El! While the results of 
Ref. (|TT|) may describe the effect of a helical potential on 
a semiconducting nanowire, it does not generally describe 
the electronic physics of a carbon nanotube. Although 
we do not study the total electronic energy, our results 
below, obtained within the long wavelength Dirac theory, 
do not show any oscillatory behavior as a function of the 
pitch. 

In Ref. dHJ), the authors studied the electronic re- 
sponse of an armchair nanotube to an applied helical 
potential. Only potentials commensurate with the nan- 
otube lattice were considered, and the high symmetry of 
the underlying armchair lattice was important in facil- 
itating calculations. This study concluded that the ex- 
ternal potential opened up small band gaps in the origi- 
nally metallic nanotube. To discover such an effect it is 
important to consider the nanotube lattice; as such, the 
continuum theory we develop below does not reproduce 
these tiny band gaps. To study a general chiral nanotube 
with an arbitrary helical potential using the method of 
Ref. (fT4| would be a formidable task, whereas the prob- 
lem is accessible within the continuum theory. The trade- 
off is that we miss higher order effects such as tiny band 
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gaps in otherwise metallic nanotubes. In general the ef- 
fects missed by the continuum theory are so small as to 
be unimportant. The corrections introduced by higher 
order considerations are discussed in the conclusion and 
in the two appendices. 

For semi-conducting tubes, we find that the band gap 
always closes under an applied, arbitrarily weak helical 
potential. For a given polymer-CNT system, the change 
in the band gap is a relatively sharply peaked function of 
pitch, with an optimum pitch maximizing the response 
to the potential. As the pitch goes to zero, the effect 
of the helical potential vanishes. As the pitch goes to 
infinity, there are two different behaviors depending on 
the strength of the applied field: for weak fields the effect 
of the chiral potential vanishes, for strong fields the band 
gap remains closed. For metallic tubes, we find that the 
Fermi velocity is a slowly decreasing function of pitch. As 
with semi-conducting tubes, as the pitch goes to zero the 
effect of helical potential vanishes. As the pitch goes to 
infinity, there are again two different behaviors depending 
on the field strength. For strong fields two small gaps, 
inversely proportional to the pitch, appear near the Fermi 
surface in addition to the gapless Fermi point. For weak 
fields no such gaps appear in the spectrum. In both cases 
the Fermi velocity is reduced by the helical potential. 

Early on it was conjectured that the underlying chi- 
rality of the CNT might determine the pitch of an ad- 
sorbed polymer^. However, it is now generally believed 
that the adsorbed species need not conform to the lat- 
tice structure of the nanotube. The structure of most 
wrapping molecules is actually incommensurate with the 
nanotube lattice, and deforming the molecule to match 
the lattice can be energetically unfavorable.^^ Simula- 
tions show that, depending on the nucleotide base se- 
quence, single stranded DNA (ssDNA) can wrap around 
a given nanotube in a left-handed or right -han ded he- 
lix, or even bond linearly along the tube.l^ Experi- 
ments show that identical strands of DNA will wrap 
with the same handedness around enantiomeric pairs 
of nanotubes.!^ Additionally, recent simulations demon- 
strate that DNA bases can adhere to a nanotube in hun- 
dreds of stable configurations. All of this is strong ev- 
idence that the underlying lattice of the nanotube plays 
little to no role in determining the structure of the CNT- 
DNA (-polymer) complex. The evidence suggests that 
the DNA (polymer) wraps in a manner determined by 
its own physical properties (chemical composition, elastic 
stiffness, size, etc.) and the CNT radius Furthermore, 
the interaction between the nanotube and the wrapping 
molecule should be adequately described by a position in- 
dependent binding energy.t^ These observations motivate 
us to ignore the discrete atomic nature of the nanotube 
and treat it as a cylinder in the continuum limit in our 
calculations below. 

The remainder of this paper is organized as follows. 
In Sec. In] we will describe our model and a convenient 
coordinate transformation. In Sec. |III| we will give ana- 
lytic and numerical solutions of the low energy spectrum 



for both metallic and semi-conducting tubes. In Sec. IV 
we use a simple model of the wrapping molecule to esti- 
mate the magnitude of the applied field and the resulting 
change in the nanotube band gap. In Sec. |V]we give a 
brief conclusion. 



II. 



MODEL SYSTEM 



As discussed in the introduction, the nanotube will be 
modeled as a continuous cylinder subjected to an exter- 
nal helical potential. By developing a long wavelength 
continuum theory we can avoid complications that arise 
when considering electrons subjected to two incommen- 
surate potentials on the scale of the lattice constant 
The low energy electronic states of the nanotube are ob- 
tained by expanding the graphene Hamiltonian around 
the K and K' points and applying appropriate boundary 
conditions. In what follows we consider the solutions 
near the K point since the response to a static helical 
potential must be the same at K and K' due to time 
reversal symmetry. 

A long wavelength continuum theory is appropriate if 
the pitch of the potential is large compared to the nan- 
otube lattice spacing. Below we will investigate the elec- 
tronic response of the nanotube as a function of the pitch 
of the potential and the nanotube circumference. The im- 
portant parameter will be the ratio of circumference to 
pitch, and we will study the response of the nanotube in 
the limits that this ratio goes to zero and as it goes to 
infinity. In the first limit we imagine fixing the nanotube 
radius and allowing the pitch to go to infinity, where our 
long wavelength description is certainly valid. In the sec- 
ond limit the long wavelength description is only valid 
if we imagine that the pitch is fixed and the radius is 
allowed to go to infinity. In the limit where the radius 
is fixed and the pitch goes to zero the underlying lattice 
structure becomes important and a continuum theory is 
no longer appropriate. 

Within the continuum theory, we note that this prob- 
lem can be reduced to a lattice periodic problem in a 
twisted coordinate system that winds with the period of 
the external potential. We choose our coordinates so that 
the X-axis lies along the tube and the y-axis wraps around 
the circumference in a counterclockwise direction (in this 
way the y coordinate agrees with the usual azimuthal an- 
gle of radial coordinates) . We let the nanotube radius be 
R and the pitch of the helix be c, as depicted in Fig. [T| 
The Dirac Hamiltonian in the presence of an external 
field becomes 



(1) 



where Vf is the Fermi velocity, a = a^x + ayj) are 
the Pauli matrices, and V{x, y) is the applied heli- 
cal potential. The envelope function is subject to 
the quasi-periodic boundary condition ^'(x,?; -t- 2ttR) = 
exp (27ri(5)^'(a;, J/), where 6 = ±1 /3 for semiconducting 
tubes and 6 = for metallic tubes .'i^Eni 
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FIG. 1: The tube on the left illustrates the geometry of our 
model system. Here, R is the radius of the nanotube and 
c is the pitch of the helical potential. On the right is an 
"unwrapped" view of the nanotube surface, with the pitch 
and the nanotube circumference, 27r7?, marked. Also on the 
right are the four relevant spatial unit vectors for this system, 
the usual i, and j/ of a Cartesian coordinate system and the 
non-orthogonal u and v used as a convenient basis in which to 
solve the Dirac equation in the presence of a helical potential. 



The helical symmetry of the external potential implies 
that V(x^y) is a function of the single variable w = 
y / R±2'iTx/ c, where the plus (minus) sign applies to left- 
handed (right-handed) hehces, and V{lo + 2tt) — V{llj). 
We will capture the essential physics of the system by in- 
vestigating the effects of the lowest Fourier mode of the 
potential. We let V{x, y) = Acos [{y — Px)/R], where A 
is the potential strength and P = ±27r i?/c is a convenient 
dimensionless parameter that characterizes the geometry 
of the helical potential, with P > for right-handed he- 
lices and P < for left-handed helices. 

In Appendix [All we investigate the effects of higher 
Fourier terms and show that for the intravalley scatter- 
ing considered below the higher Fourier terms are unim- 
portant. However, higher Fourier terms can vary rapidly 
on the scale of a lattice constant, and rapidly varying 
potentials can lead to significant intervalley scattering.ESl 
Intervalley scattering involves a large momentum trans- 
fer and in general the Coulomb potential is unscreened 
at large momenta. Such an unscreened interaction gives 
higher order Fourier coefficients that are larger than the 
lowest screened Fourier coefficient, and we might expect 
that intervalley scattering will be at least as important 
as intravalley scattering. However, in Appendix [X]2 we 
show that in general the intervalley scattering matrix el- 
ements are either kinematically forbidden or small com- 
pared with the first screened coefficient. In the very few 
cases where intervalley scattering may be important, we 
show that the matrix elements are about the same size 
as the intravalley elements and do not appreciably alter 
any of the results below. 

The isotropy of the Dirac equation allows us to choose 
any two unit vectors as our basis vectors in the tangent 
plane of the tube. Here, it is convenient to take as a basis 
two unit vectors, u and v, such that u points along an 
equipotential and v follows the usual azimuthal coordi- 



nate. Referring to Fig.jT] we see that 

X Py 



V ^ y. 



(2) 



With this choice of basis vectors a general vector in the 
plane is written as r = uu -I- vv. The {u^v) coordinates 
may be obtained from the Cartesian (x, y) coordinates 

by 



P2 
-P 1 



V 



(3) 



In this basis the helical potential is a function of only the 
V coordinate, V{v) = Acos (v/R), and the Dirac equation 
becomes 



- ihvf (^Vl + P^(Jxd^, + (ay - Pa^)dyj * 

-fAcos(u/P)«' = (4) 

With this choice of spatial basis vectors we are forced 
to choose non-orthogonal reciprocal lattice basis vectors. 
We choose our reciprocal basis vectors, (f„ and such 
that for a general wavevector, k — kuqu + k^q^, the dot 
product with a general spatial vector is given by fc • r = 
kuU + kyV. This is accomplis hed usi ng the reciprocal 
space (non-unit) vectors Qu = Vl + P^x and (jv = —Px+ 
y. With this basis the {ku,ky) coordinates are obtained 
from the usual (fcj,, ky) coordinates by the transformation 



1 



1 



P 



(5) 



If we now write ^'(u, u) = exjp (ikyujip^v), the Dirac 
equation becomes 



1 + P^ik^R)a, - iR{ay - Pa,)dyj ^ 

+a cos {v / R)^p — eip, (6) 



where a = A/Aq and e = E/Aq are dimensionless mea- 
sures of the potential strength and the energy, respec- 
tively, and Aq = hvf/R is a convenient unit of en- 
ergy. The above Hamiltonian has the property that 
Hjii—ky, —P) — (TyHxiky^, P)ay, which shows that the 
spectrum satisfies e(A:„, P) = e{—ku, —P)- We may there- 
fore only consider right-handed helices and take P > 
in all that follows. 

When P = Eq. (|6| also describes a nanotube im- 
mersed in a constant perpendicular electric field, a sys- 
tem studied previously by Novikov and Levitov.'SllIl 
They found that for semi-conducting tubes the band gap 
was unaffected by the applied field until a critical field 
strength was reached, above which the band gap closes. 
For metallic tubes there exists a critical field above which 
the Fermi velocity changes sign and the Fermi surface 
fractures. We will investigate similar effects for P 7^ 
below. 
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Novikov and Levitov used a chiral gauge transforma- 
tion to show that, for the P = system, the spectrum at 
/cj; = is unaffected by the appHed field. For the P ^ 
system we employ a similar transformation. 



T = exp ( ia{(Ty — P<Tx) 



sin (v/R) 



(7) 



which converts the Hamiltonian to 

H'j, = THkT-' 

- -iR{ay - Pa,)d, + + P^{KR)T(7^T-\ 

(8) 

Thus, the spectrum of Eq. (6]) at fc^, = is unaffected by 
the applied helical potential. 

When fc„ = solutions of Eq. ^ are of the form 



^{v) = f{v) 



where /(w) is independent of a. This shows that the 
applied field only alters the phase of the wavefunction 
at = 0, and we use this fact to elucidate another 
protected quantity in the spectrum of Eq. ^ . The k^- 
space energy gradient is given by 



dku 



(10) 



provided both de/dku and d\ip)/dku are well defined. 
Note that in general these derivatives are not well defined 
at degenerate points in the spectrum. For non-degenerate 
fc„ points, if we write the wavefunction as 



then Eq. ( 10 1 becomes 

^ = R^l + / dv {(t>*e + 



(11) 



(12) 



By Eq. ([9]), the right hand side is independent of a at 
ku = 0, which shows that de/ dku\ku=o is unaffected by 
the external helical potential. We therefore find that the 
longitudinal velocity, given by 



1 dE 
h dkx 



Ao 



de 



hVl + P2 dku ' 
is unaffected by the helical potential at fc„ — 0. 



(13) 



III. 



SOLUTIONS 



Zero Field Solutions 



The helical band structure of Eq. ^ is different than 
that produced by the ordinary massless Dirac theory on 




FIG. 2: The three highest valence bands and the three lowest 
conduction bands of a semiconducting nanotube {5 = 1/3) 
plotted as a function of kuR- The energy bands are plotted 
for a free nanotube (a — 0) studied in a twisted coordinate 
system appropriate for an apphed potential with P = 1. The 
allowed energies, like all measurable quantities, are unaffected 
by the change in coordinate system. The apparent differences 
between this spectrum and the usual free nanotube spectrum, 
such as the shift in the location of band minima, arise because 
the energies are plotted as a function of kuR, defined for the 
twisted coordinate system, and not the usual k^R. 



the cylinder both because of our choice of nonorthogonal 
coordinates and because of the applied helical potential. 
Only the latter change is physical, and in order to sep- 
arate the two effects it is necessary to examine the field 
free solutions of Eq. ([6| . The appropriate boundary con- 
dition on ip{v) is obtained by noting that v = y, so that 
ip{v + 2ttR) — exp {2Trid)tlj{v). The field free energies are 
given by 



Vl + P2(fc„i?) - P{m + S) +(to + <5)2, 

(14) 

where the band index m is any integer and a = ±1 
denotes a conduction or valence band. The first few 
energy levels for a semiconducting CNT are plotted in 
Fig. |2] The band gaps, like all observable quantities, 
are unchanged by our choice of coordinate system. How- 
ever, there are two noticeable differences when the free 
nanotube spectrum is plotted as a function of fc„i?, de- 
fined for the twisted coordinated system, compared to 
the spectrum when plotted as a function of the usual 
kxR. First, the band minima are shifted to positive 
(negative) fc„ values for m > (to < 0). The loca- 
tions of the new band minima are given by {kuR)min — 
P{m + S) + Secondly, as k^R oo, the slope of 
each band goes as \/l + P'^ . 

An intriguing feature of the band structure of semi- 
conducting CNTs is that the bands are shifted so that 
the magnitude of the slope of all bands is the same at 



0, namely 



PR. As these 
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FIG. 3: The lowest conduction band and highest valence 
band as function of feu 7? for a semi-conducting CNT with 
(5 = 1/3 and P = 0.1, drawn for several values of a as labeled 
in the figure. The inset is an expanded view of the lowest 
conduction bands in the area around fc^ii = 0, and shows 
that the band gap closes for every non-zero a considered here. 



are the protected locations in the spectrum, it is those 
states with A:„ = and with longitudinal velocity v = 
±Pw//vTT~P^ that are unaffected by the applied po- 
tential. The same holds true for the energy bands of 
metallic tubes except for the m = bands. These bands, 
with energy tafi — ol\/\ H- P'^{kuR), remain degenerate 
at ku = 0- From the analysis at the end of Sec. [iT] the 
slope of these bands may be changed by the applied field. 
These changes will be investigated below. 



B. Non-zero Field Solutions 



FIG. 4: The lowest conduction band and highest valence 
band as function of kuR for a semi-conducting CNT with 
S — 1/3 and P = 1.0, drawn for several values of a as labeled 
in the figure. 




1. Semiconducting Nanotubes 

The spectrum of Eq. ^ for semi-conducting nan- 
otubes was determined numerically for many values of 
the field strength, a, and the dimensionless geometric pa- 
rameter, P. A representative sample of results is shown 
in Figs. [3]- [5] where we plot the lowest conduction band 
and highest valence band for three different values of P 
and three different non-zero values of a. The most strik- 
ing result of these solutions is that, unlike the P — 
case, the band gap closes for any value of a. This is 
shown in the inset in Fig. [3] and will be confirmed an- 
alytically below. Although there is no critical field to 
be applied before the band gap closes, our results sug- 
gest that, at least for P < 0.2, there exists a critical field 
above which a second local minimum appears in the band 
structure. This is most clearly seen in Fig. |3] where the 
a — 1.0 band develops a second local minimum at about 
kuR~ —0.48. For larger P there is significant mixing of 
the m = and m = — 1 bands in the region where a sec- 
ond local minimum would develop, and the existence of 



FIG. 5: The lowest conduction band and highest valence 
band as function of kuR for a semi-conducting CNT with 
5=1/3 and P = 2.0, drawn for several values of a as labeled 
in the figure. 



a critical field in these cases is less certain. Nevertheless, 
for all values of a and P > there exists only one global 
minimum. The location of the band minimum, (A;u^)mini 
increases as a increases. 

The fractional change in the band gap is given by 
F = [A(P,a) - A(P,0)]/A(P,0), where A(P, a) is the 
minimum value of the conduction band for a given P 
and a, and A(P, 0) = 1/3 is the minimum value of the 
unperturbed conduction band. F is plotted in Fig. [6] as 
a function of P for several values of a. For all a, F has a 
maximum for a P on the order of unity and decays rela- 
tively rapidly to zero as P — *■ 00. This may also be seen 
by comparing Figs. |4] and [S] which show that doubling 
P significantly decreases the effect of the applied field. 
This behavior is most easily understood by considering 
the tight-binding model. The limit P — > cx) corresponds 
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FIG. 6: The fractional change in the band gap, F, plotted 
as a function of P for several values of a. The band gap is 
unchanged as P ^ oo, as discussed in the text. As P ^ 
there are two behaviors: the band gap is unchanged if a < Sc, 
but closes for a > Uc- For all a there exists a P on the order 
of unity that produces the largest change in the band gap. 



to tciking the pitch to zero, at which point the helix col- 
lapses into a uniform cylinder. The potential is no longer 
spatially varying on the tube surface, and its only effect 
is to provide a uniform background potential. Such a 
uniform change may be eliminated by redefining the zero 
of energy, and thus cannot affect the band gap. 

As P ^ the system asymptotically approaches the 
system studied by Novikov and Levitov, and their results 
apply. At P = the band gap is unaffected by the ap- 
plied field if a < Oj, « 0.6215, and it closes if a > Uc- 
These two behaviors are evident in Fig.|6j for a < Qc, F 
rapidly approaches zero as P ^ 0, while for a > ac, F 
approaches an a-dependent constant less than zero. 

To investigate the size of the band gap for small a we 
employ non-degenerate perturbation theory. The per- 
turbation connects band (a,m) to bands {zt,m+ 1) and 
(±,m — 1), and the perturbed energies are given by 



(i 



2 G 



(15) 



where e° is given by Eq. ( 14 1 and 



A[P{kuR)-Vm^{m + S)] -(1 + P2) 



(16) 

Because of the many degeneracies in the spectrum at 
high energy (see Fig. |2]), the effects of the chiral poten- 
tial on the entire spectrum cannot be studied within a 
non-degenerate perturbation theory. However, there is 
never a degeneracy at the band minimum of the m = 
band, and Eq. (15) accurately reproduces the band gaps 



found in numerical results. The fractional change in the 
band gap may be estimated by evaluating Eq. ( 15 1 at the 



Fa,o = a 



(I) 



2 4P2|5| 

4,52- (1 + P2)2- 



(17) 



Notice that as (1 + P^) > 2\6\ for all P, the band gap 
always closes. The actual band gap is in fact slightly 
smaller than predicted by Eq. (17 1 because the unper- 



turbed band minimum is not the true band minimum. 
Nevertheless, the conclusion that the band gap closes for 
any a > remains valid. 

To investigate the existence of a critical field and the 
development of a second local minimum we employ non- 
degenerate perturbation theory in small fc„P. The per- 
turbation connects conduction bands to valence bands, 
and the perturbed energies are given by 



aim - 



(5|\/l-t-P2 - a sgn(m + (5)P(fc„P) 



+asgri(m + S) , . 

Vl + P2 ■ 

where sgn(x) is the sign function and 



H„ 



00 

E 

n— — 00 



r 



2a 



2{m + S) 



Vl + PV 4(m-h(5)2 



(18) 



(19) 



1/3, Hq switches sign from positive to negative 
For a 



For S 

at flc w 0.6215\/1 + P^- For a > ttc the curvature at 
ku = is negative and a second local minimum develops 
to the left of the origin, as deduced from the numerical 
results. 



2. Metallic Tubes 

The low energy spectrum of Eq. ^ for metallic tubes 
is shown in Figs. |7] and [8] for two different values of a 
and several values of P. As discovered by Novikov and 
Levitov, for P = 0, as the field is turned on the Fermi 
velocity decreases and the band acquires a non-zero cur- 
vature away from fc„ = 0. The Fermi velocity goes to zero 
at a critical field strength, Ccm ~ 1.203. When a > acm 
the Fermi velocity switches sign and the Fermi surface 
fractures, as shown in Fig. |8] 

The spectrum for fixed a < is plotted in Fig. [7] 
for several P values. As P increases the slope of the 
energy bands at fc„ = increases. However, note that 
as Vf — de/dkx — (l/Vl + P^)de/dku, the increase in 
the slope does not necessarily translate into an increase 
in the Fermi velocity. In fact, for very large P the slope 
increases linearly in P, which implies an unchanged Fermi 
velocity in the limit P ^ 00. This effect could have 
been anticipated from the analysis of the semiconducting 
band gap, where we argued that the effect of the chiral 
potential must disappear in the limit P —i- 00. As P 
increases the higher energy states mix with the m = ±1 
bands and two symmetric local minima develop in the 
lowest band. 



7 



1.0- 



0.5- 



An 



-0.5- 



-1.0- 



P=1.0 



_ P=2.0 



,.P=0.5 



P=0.25 ■ 



-1.0 



-0.5 







0.5 



1.0 



FIG. 7: The low energy bands of a metallic nanotube as a 
function of kuR for a = 1.0 < Ocm for several values of P 
as labeled in the figure. The primary effect of the field is to 
modify the Fermi velocity. 




FIG. 8: The low energy bands of a metallic nanotube as a 
function of kuR for a = 1.5 > acm for several values of P as 
labeled in the figure. For P = the Fermi surface fractures, 
for P > the Fermi surface is unchanged. 



The spectrum for fixed a > acm is plotted in Fig. [8] 
where it is seen that the effect of a non-zero P is to open 
up gaps where the P = bands crossed the Fermi surface 
at non-zero fc„ . For any P > the Fermi surface does not 
fracture, and for this property P = is a singular limit. 
As P increases the slope at A:„ = increases and the local 
minima at the induced band gaps at first flatten, then 
reappear as the lowest band begins to mix with the m — 
±1 bands. As P — *■ oo the effects of the chiral potential 
disappear and the a > acm system maps directly onto 
the a < acm system. 

We study the change in the Fermi velocity using de- 
generate perturbation theory on the two states at A:„ = 0. 
The perturbed energies are given by 



P2 



2a/Vl + p2 (/c„P) 



(20) 



For P = the Fermi velocity vanishes when 2a = /ii « 
2.405, the first zero of the Bessel function. When P > 
the Fermi velocity is strictly non-zero. The slope of the 
bands is an oscillatory function of a, but as Jo{x) < 1, 
this oscillatory behavior will be small for large P. As 
P oo the slope goes as P, in agreement with the nu- 
merical results above. 



IV. 



THE FIELD STRENGTH 



To a first approximation, the DNA and polymers used 
to wrap nanotubes may be considered infinitely long, 
helically wrapped line charges. In order to estimate 
the magnitude of the field strength we model the wrap- 
ping molecule as a helical ribbon of width Z, radius Pi, 
and pitch c, with surface charge density ctq (see Fig. |9|. 
The field strength is obtained by evaluating the leading 
Fourier coefficient of this charge configuration on the nan- 
otube surface, r — R. At the end of this calculation we 
recover the line charge model by taking the limit I 
with the linear charge density, A = Iuq, held fixed. In 
this section we work in SI units. 

In a long wavelength theory of the applied fields, we 
break space into three macroscopic regions characterized 
by different dielectric constants, as shown in Fig. |9] Re- 
gion 1, r > Pi, is composed of material outside the heli- 
cal line charge, the aqueous solution or other solvent, or 
possibly vacuum. Region 2, Ri > r > R, is composed of 
the material between the line charge and the nanotube, 
including DNA bases or uncharged regions of a wrapping 
polymer. Region 3, P > r > 0, is composed of the nan- 
otube itself. Denote the dielectric constant in the zth 
region by e^. 

The helical symmetry of the charge distribution im- 
plies that the potential can be expanded in cylindrical 
coordinates as 

oo 

F,(r, </>, z) = In (r/r^^^) + ^ (AW/„(27rnr/c) 

ri=l 

+ BlpKn{2Trnr/c)'^ cos[n(0 - 27rz/c)], 

(21) 

where i — 1,2,3 denotes the three regions of space and 



Al:\Bi'\ and Tq are expansion coefficients to be deter- 
mined. The coefficients are determined using the usual 
electrostatic boundary conditions,!^ 



9V> 



dr 



V>{r = r^) = V<(r = r^) 



dr 



a{r^n)/eo, (22) 



where the boundary is located at r = r.^, V> and e> 
(V^ and e<) denote the potential and dielectric constant 
for r > ri (r < r^), and a is the charge density at the 
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FIG. 9: On the left is a depiction of our model, a tube of 
radius 7? surrounded by a helical ribbon of radius Ri, width I, 
and surface charge density uq . On the right is a cross-sectional 
view of the system. The numbers refer to the three regions of 
space with different dielectric constants, as discussed in the 
text. 



boundary. The charge density at r — Ri is given by 
a{(j), z) — ° ° + — sin(n7ra;o/c) cos[n(0 — 27rz/c)], 

n—l 

(23) 

where xq = ly^l + Pi/ Pi and Pi = 2'kRi/c. 

The field strength is obtained from the Fourier expan- 
sion by a 



eA^^''/i(P)/Ao, where the electron charge 
is — e. After taking the limit as the ribbon goes to a line 
charge, we find 



eA 



-W{Ri,R) = CW{Ri,R), 



where C is independent of R and i?i , and 

'63 \ RU{Pi,P) 



(24) 



with 



S{x,y) 



eiy K[iPi)SiPi,P) - {e2/ei)Ki{Pi)T{Pi,Py 

(25) 



Vl -1- a;2 


Ki{x)Ii{v) 


(Kiiv) 


K'M\ 


X 


1 - (£3/62) 


\ h{y) 


I'M ) 


Ki{x) + 


h{x) 


^3 Ki{y) 


K[{y) 


1 - (£3/62) 


\<^2 hiv) 


I'liy) 


K[{x) + 


I[{x) 


^3 Ki{y) 


K'M 


1 - (£3/62) 


li{y) 


I'M 



(26) 

To determine C we consider specifically the case of 
an ssDNA-CNT complex. The linear charge density of 
ssDNA is obtained by assuming each phosphate group 
on the backbone carries a chaxge of — e, which gives 
A ~ — 1.5e/nm. The dielectric constant of a CNT is 
obtained by using a result from Refs. (j24| and (l2T|l . that 
when immersed in a perpendicular electric field the ra- 
tio of the field strength inside the tube to the applied 




20 30 
c (nm) 



FIG. 10: The function W{R, c, ei) as a function of c for two 
values of R and two values of e\. 



field is given by E/Eq = 1/5, independent of R. If the 
nanotube is modeled as a uniform solid cylinder with di- 
electric constant 63 then E/Eq = 2/(e3 + 1), which gives 
£3 — 9. Using these results and « 8 x 10^ m/s, we find 
C « -1.8 nm-i. 

The distance between the ssDNA backbone and CNT 
surface is independent of nanotube radius, and is given 
hy Ri — R ~ 0.6 nm. The dielectric constant for region 
2 is difficult to approximate as it is composed of ran- 
dom DNA bases, solvent molecules and other dissolved 
species. However, provided the pitch is large compared 
to the nanotube radius the dominant screening effects 
should come from regions 1 and 3. We therefore ignore 
screening in region 2 and set 62 = 1. 

With these approximations W = W{R,c,ei). The 
asymptotic behavior of W for R ^ I nm and ei > 5 
is ~ woR/ei, where wq is a function of c. The be- 
havior of ly as a function of c is shown in Fig. [TO] for 
two values of R and ei. W goes to zero as c ^ 0, and 
goes asymptotically to a constant value as c — > 00. Phys- 
ically realizable values of c are in the 1 — 60 nm range,'^ 
and in this region W is monotonically decreasing but its 
behavior is not given by a simple scaling relation. 

To obtain an estimate of the field strength we use a 
CNT with R = 0.5 nm and a solution with ei = 80, 
which is approximately the dielectric constant of water. 
In Table [l] we calculate the field strength and the corre- 
sponding fractional change in the band gap for several 
values of the pitch. As c increases a asymptotically ap- 
proaches 1.1 X 10~^, while P monotonically decreases to 
zero. From Fig. [6] we expect the magnitude of F will 
be peaked for a P on the order of unity and decay to 
zero for very small and very large pitches. This behav- 
ior is confirmed in Table [T] where c — 5 nm produces 
the largest change in the band gap. Notice that as the 
pitch increases the field strength continues to increase, 
but the corresponding decrease in P reduces the effect of 
the applied potential. 

For the largest fractional change found here, the dif- 
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TABLE I: Calculated values of P, the field strength, a, and 
the fractional change in the band gap, F, for various values of 
the pitch, c. The other parameters of the system are described 
in the text. 



c (nm) 


P 


a 


F 


1 


3.14 


6.6 X 10"" 


-3.6 X lO"** 


2 


1.57 


3.7 X 10"=* 


-2.9 X lO"** 


5 


0.63 


7.6 X 10"^ 


-1.5 X 10"^ 


10 


0.31 


8.8 X 10"^ 


-1.0 X 10"^ 


25 


0.13 


1.0 X 10"^ 


-2.7 X 10~^ 


50 


0.063 


1.0 X 10"^ 


-7.7 X 10"^ 



ference between the original and perturbed band gap is 
about 0.01 meV. We do not expect a change of this scale 
to be readily measurable in transport or optical experi- 
ments. There are several possible ways one might con- 
sider to increase the size of this effect. The first is to use 
large radius nanotubes, which will increase the magni- 
tude of the applied chiral potential. This approach has 
two problems. Fundamentally, as R is increased P is 
also increased, and the size of the effect decreases with 
increasing P. Practically, single-walled nanotubes are 
generally no larger than 1.0 nm in diameter, with larger 
radius tubes unstable to collapse.'^ The second method 
is to tune c so that for a given radius nanotube, P is 
slightly less than unity where the effect of the helical po- 
tential is largest. For shorter DNA strands a systematic 
study of pitch with varying bases is possible,^ but for a 
long DNA strand with hundreds or thousands of bases 
such a study is impractical. A third method is to change 
the wrapping polymer to one with a larger linear charge 
density or to one that rests closer to the nanotube sur- 
face. The linear charge density of DNA is already rather 
large and simply doubling or tripling the linear charge 
density will not increase the size of the effect enough for 
easy observation. The field strength increases exponen- 
tially as the charge approaches the nanotube surface, but 
the decay length is large relative to typical molecular dis- 
tances. A fourth method is to change the environment 
of the DNA-CNT complex. The dielectric constant of 
water is unusually large and essentially reduces the field 
strength by a factor of 80. Most organic solvents have di- 
electric constants in the range 1 — 5, which would increase 
the field strengths calculated above by factors of 16 — 80. 
The DNA-CNT complex is remarkably stable, and it may 
be possible to create the complex in aqueous solution and 
then transfer it to an organic solvent. Other wrapping 
molecules hav e alre ady been used to dissolve CNTs in or- 
ganic solvents ,1^^^ and one of these might produce field 
strengths strong enough to produce a measurable effect. 



V. CONCLUSION 

We have investigated the effect of a helical potential 
on semi-conducting and metallic nanotubes. For semi- 



conducting nanotubes the band gap closes for any non- 
zero field strength. The size of the effect is determined by 
both the field strength, a, and a dimensionless geomet- 
rical factor P, which is the ratio of the circumference of 
the nanotube to the pitch of the helix. For each a there 
exists an optimal P that produces the biggest change in 
the band gap. For metallic tubes, the helical potential 
decreases the Fermi velocity but does not fracture the 
Fermi surface. Under typical conditions the effect of the 
helical potential is probably unobservable, but we un- 
derstand the scaling relationship between the size of the 
band gap and every control parameter, so in theory it 
is possible to design a system where these effects would 
be detectable in an optical experiment. In some ways 
the small size of the effect is encouraging, as it means 
that measurements on helically wrapped CNTs give re- 
sults that are nearly identical to those of pristine CNTs. 
Also, a helically wrapped CNT can be substituted for a 
pristine CNT in almost any application, which may make 
device construction easier. 

The theory developed here ignores higher order correc- 
tions, such as curvature effects and higher Fourier terms 
in the potential. Simply expanding the tight-binding 
graphene Hamiltonian to the next lowest gradient order 
generates terms that break the chiral gauge symmetry 
and introduce a correction to the spectrum at fc„ = 0. 
Such corrections will modify the functional form of our 
analytic expressions; for example, curvature effects will 
introduce a chiral angle dependent band gap. However, 
curvature effects and higher order expansion terms in- 
troduce corrections that are small compared to the ener- 
gies of interest. The theory developed here also ignores 
exciton effects, which are known to be large in CNTs. 
Exciton effects in optical experiments may alter the nu- 
merical values obtained here by 20—30%, but they should 
not significantly change our general results or alter the 
conclusions outlined in the previous paragraph. 

This work was supported by the Department of Energy 
under grant DE-FG-84ER45118. 



APPENDIX A: HIGHER FOURIER 
COEFFICIENTS 

A general chiral potential on the nanotube may be 
written as 

oo 

V{f) = Y,Vr.COs{Qn-f), (Al) 

n=l 

where Q„ = ±27rn/c x+n/R y. In our calculations above 
we only retained the first Fourier term in order to cap- 
ture the essential physics of the system. In this Appendix 
we investigate the effect of including the higher Fourier 
terms in the Hamiltonian. Our perturbation calcula- 
tions above only included contributions from intravalley 
scattering matrix elements, specifically, scattering from a 
state near the K point to another state near the K point. 
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TABLE II: Calculated values of higher field strength Fourier 
coefficients for the nanotube system studied in Sec. |IV| 



c (nm) 


ai 








a/1,1 a\ 




2 


0.0037 


0.063 


5.4 X 10" 


-•i 


5.2 X 10" 


-4 


5 


0.0075 


0.16 


3.7 X 10" 


^2 


9.5 X 10" 


'3 


10 


0.0088 


0.20 


5.6 X 10" 


-2 


1.8 X 10" 


-2 


50 


0.0104 


0.20 


5.9 X 10" 


^2 


2.0 X 10" 


-2 



In section [Xjl we show that in all cases higher Fourier 
coefficients add a negligible correction to the intravalley 
matrix elements. 

The higher Fourier terms vary rapidly on the scale of 
a lattice constant and contribute to scattering with a 
large momentum transfer. Such scattering may connect 
states at K with those at K', so we must also consider in- 
tervalley scattering matrix elements in our perturbation 
expansion. At large momentum transfer the Coulomb 
interaction is unscreened and the Fourier coefficients in 



Eq. ( Al I will be significantly larger than the correspond- 
ing screened coefficients. Nevertheless, we show in sec- 
tion |Aj2 that in almost all cases the intervalley scatter- 
ing matrix elements can be neglected compared to the 
intravalley scattering matrix elements. In a few systems 
with small radius nanotubes, the intervalley matrix ele- 
ments may be nearly as large as the intravalley scattering 
matrix elements. We discuss these few cases below. 



where a„ = — el^/Ao. The coefficient a\ was evaluated 
in Sec. IIVI The nth Fourier coefficient is obtained in the 
same manner and is given by equations very similar to 



Eqs. (24 - 26 1, but involving ln(nx) and Kn(nx) instead 
of and K\{x). If this full potential is used to calcu- 
late the intravalley matrix elements, then the nth Fourier 
term connects the m = band to the m = ±n bands. 
When squared and summed to obtain the second order 
perturbation to the energy, all of the cross terms vanish 
and we are left with 



cia\ 



c-ial 



cial [I + {ailaxfl2 + {a^i/aif/i + 



(A3) 



where the Ci are constants (independent of the field 
strength) and we have approximated c„ « ci/n, the fac- 
tor of 1 /n coming from the energy denominator. In Ta- 
ble |lT] we Ust a^/ai, for i = 2 - 4, for the DNA-CNT 



system studied in Sec. IV Using 02/01 < 0.2, we find 
that including the second Fourier term only changes the 
second order energy shift by about 2%. The energy shift 
was already much less than the unperturbed energy, and 
a 2% change to such a small shift can obviously be ne- 
glected. 



1. Intravalley Scattering 

The dimensionless potential energy appearing in the 2. Intervalley Scattering 
Hamiltonian is given by 

00 The effects of intervalley scattering can be included 

U{r) — a„cos((3„ • r), (A2) with an effective Hamiltonian derived in Ref. ()20| . 



-ihvf {dx + idy) 
-cj*e-*^u^(r) 



v 







-ihvf {dx - idy) 
UB{r) 




J6~, 



UB{r) 





UA{r) 
-ihvf (dx - idy) 





e-^'uUr) 
-ihvf {dx + idy) 
UB{r) 



(A4) 



where the states are written in the KA, KB, K' A, K' B effective potentials in Eq. (A4| are given by 
basis, 9 is the chiral angle, and uj = exp(27ri/3). The 



UA{r) 


Ra 


Ra)U{Ra). 


UB{r) 


Rb 


Rb)U{Rb), 


UA{r) 


Ra 


RA)e^(^~^y^-U{RA), 


UB{r) 


Rb 


RB)e'^^-^'>^^U{RB), (A5) 
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TABLE III: Calculated ratios of the nth unscreened Fourier coefficient to the first screened Fourier coefficient for the nanotube 
system studied in Sec. |IV| 



c (nm) 






4"V4^^ 


4"V4^' 


4"V4^' 


4"V4^' 


a("V4^' 


2 


190 


13 


1.1 


0.10 


1.1 X 10"^ 


1.2 X lO""* 


1.3 X IQ-^ 


5 


200 


35 


7.8 


2.0 


0.52 


0.15 


4.2 X 10"^ 


10 


190 


41 


12 


3.6 


1.2 


0.42 


0.15 


50 


170 


40 


12 


4.0 


1.5 


0.56 


0.22 



27r 

3a 
27r 



where C/(r) is given by Eq. (A2), Rj^ = R + {Rb = angle and are given by 
R + Tb) denotes a site on the A (B) sublattice, and 
g{r — R) is a function peaked at R with a width of about 
a lattice constant and normalized so that J^Rdi^) ^ 1- 
The effective potentials contain information about the 
external chiral potential and the underlying graphene 
lattice, and in general do not share the same symme- 
try as the applied potential. Furthermore, Ui and Ui do 
not share the same symmetry, and there is no coordinate 
transformation analogous to Eqs. ^ and ([3| to express 
these potentials in terms of a single variable. 



K[ = ^ |(^cos6' + \/3sin6'^ i+ (^\/3cos6'-sin6l^ y| 



K2 = ^|(cos6l-\/3sine')x- (^sin^ + A/Scose^yj 

^ 47r 47r 

= -— cos6'x-|- — sin6'y. (A8) 
3a 3a ^ ^ 



To remain in the low energy regime we require |p| < 
1/a, which implies that the only terms contributing to 
the sum in Eq. ( A7 ) will have 



Second order perturbation theory shows that the off- 
diagonal matrix elements contribute in the same way as 
the diagonal matrix elements (that is, we have to com- 
pute \{K\u\K')\^/{Ek - Ek'))- Thus, to get an esti- 
mate for the energy shift due to the off-diagonal terms 
we need only calculate the largest Fourier coefficient of 
those terms. As mentioned in Sec. |TT] the Coulomb po- 
tential is nearly unscreened at high momentum transfer, 
and the unscreened Fourier coefficients of the potential 
can be much larger than the screened coefficients. In Ta- 
ble 



III 



we list al"-'/a^' 



the ratio of the nth unscreened 
Fourier coefficient to the first screened Fourier coefficient, 
for the nanotube system studied in Sec. |IV[ From the size 
of the coefficients it appears intervalley scattering could 
be more important than intravalley scattering. 

The off-diagonal matrix elements are evaluated as 



{K,k^,m\u^\K',k'^,n) oc X! X! '^P-Q„+i?-i?',G' 

715^0 G 

(A6) 

where p = [kx — k'^)x + (to — n -\- 26) /R y, the sum on 
n is over all integers except 0, and the sum on G is over 
all reciprocal lattice vectors of the graphene lattice. This 
may be rewritten as 



{K, kx,m\u^\K',k'^,n) cx X! '^P-Q„:i?," 

Ti^O K'- 



(A7) 



where the sum is over all K' points in reciprocal space. 
Because of the magnitudes of p and Qn, only the K' 
points on the edge of the first Brillouin may contribute 
to this sum. These points are a function of the chiral 



\Qn~Kl\<l/a. 



(A9) 



The expression on the left-hand side of Eq. ( A9 1 depends 
on n, c, R, 0, and on the particular choice of K' point. 
We undertook a systematic search of the relevant pa- 
rameter space to determine when Eq. (A9l was satis- 
fied. The search was limited to —6 < n < 6 because 



it is clear from Table III that the coefficients of higher 
Fourier terms will always be negligible. Preliminary in- 
vestigations showed it was sufficient to restrict the chiral 
angle to 6 — 0,±7r/12 or ±7r/6. For each of these 180 
cases we plotted |Q„ — K^\ as a function of c and R to 
find where the inequality (A9l is satisfied. The search 



was restricted to physical values of the pitch and the ra- 
dius, 1 nm < c < 60 nm and 0.3 nm < R < 1.0 nm. We 
found 26 cases where there was any region in the (c, R) 
plane where the inequality was satisfied. 

For each of these regions we evaluated the ratio 
ari^ /oi^ to determine if the size of the intervalley scat- 
tering matrix element is comparable to the size of the 
intravalley scattering matrix element. As a conservative 
estimate the ratio was considered significant if 



la\"> > 0.75. 



(AlO) 



There are 7 cases where the intervalley matrix element 



is kinematically allowed (inequality ( A9 ) is satisfied) and 



the unscreened coefficient is large enough to be consid- 
ered significant. These 7 cases are listed in Table |IV| 
where it is seen that both inequalities are satisfied only 
for very small radius nanotubes, and then only in a nar- 
row range of radii. Likewise, the pitch must generally 
be restricted to a narrow range. In an arbitrary sam- 
ple of nanotubes the fraction of tubes that satisfy such 
restrictions will be small, and the intravalley scattering 
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TABLE IV; This table shows the 7 cases where both inequahties (A9l and (AlOl are satisfied. The range of c and R over 



which the inequahties are satisfied is given, along with the ratio of the unscreened coefficient to the screened coefficient and 
the fractional change in the energy shift. In the first column i refers to a specific Kl point. 



Range of c (nm) 



Range of R (nm) 



ial/aiy/\n\ 



3 
4 
4 
5 
5 

-4 

-5 



n/12 
-7r/12 


-7v/6 
-n/12 
7r/6 
7r/6 



1.9 - 2.2 

3- 13 

4- 5.5 
unrestricted 

15-20 

> 20 

> 15 



0.30-0.35 
0.30-0.32 
0.30 - 0.33 
0.37-0.38 
0.37-0.38 
0.30-0.31 
0.37-0.38 



0.75- 1.2 
0.75 - 1.75 
0.75- 1.3 
0.75 - 0.90 
0.75 - 0.80 

~ 1.75 
0.75 - 0.80 



0.19-0.48 
0.14-0.76 
0.14-0.42 
0.11 - 0.16 
0.11 - 0.13 

~ 0.76 
0.11 - 0.13 



effects will dominate the intervalley scattering effects in 
any measurement. 

It may be possible to prepare a sample in one of the 
allowed regions in Table |IV[ but even then intervalley 
effects will be no more important than the intravalley 
scattering already calculated. The energy shifts due to 
intravalley scattering are so small compared to the un- 
perturbed energy that adding an additional shift of the 
same order of magnitude will not make the effects of a 



chiral potential any easier to observe. Furthermore, the 
additional shifts due to intervalley scattering are small 
enough that they do not alter any of the conclusions from 
the main body of the paper. Lastly, note that these re- 
sults were calculated assuming no screening. If there is 
any residual screening either by the nanotube or by the 
environment, then the shifts due to intervalley scattering 
will be further suppressed. 
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